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Abstract 

We present release 1.3 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an international open source 
software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models 
such as quantum magnets, lattice bosons, and strongly correlated fermion systems. Development is centered on common XML 
and binary data formats, on libraries to simplify and speed up code development, and on full-featured simulation programs. The 
programs enable non-experts to start carrying out numerical simulations by providing basic implementations of the important 
algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble 
simulations, exact and full diagonalization (ED), as well as the density matrix renormalization group (DMRG). Changes in the 
new release include a DMRG program for interacting models, support for translation symmetries in the diagonalization programs, 
the ability to define custom measurement operators, and support for inhomogeneous systems, such as lattice models with traps. 
The software is available from our web server at http://alps.comp-phys.org/. 
© 2008 Elsevier B.V. All rights reserved. 

PACS: 02.70.-c; 75.40.Mg 
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1. Introduction 

In this paper we present release 1.3 of the ALPS project 
(Algorithms and Libraries for Physics Simulations), an 
open source software development project for strongly cor- 
related lattice models. We will present a short overview and 
focus on new features compared to the previous releases 
[1,2]. 

Quantum fluctuations and competing interactions in 
quantum many body systems lead to unusual and exciting 
properties of strongly correlated materials such as quan- 
tum magnetism [3], high temperature superconductivity 
[4], heavy fermion behavior [5], and topological quantum 
order [6]. The same strong interactions make accurate an- 
alytical treatments hard and direct numerical simulations 
are essential to increase our understanding of the unusual 
properties of these systems. 

The last decade has seen tremendous progress in the de- 
velopment of algorithms. Speedups of many orders of mag- 
nitude have been achieved [7,8,9,10,11,12]. These advances 
come at the cost of substantially increased algorithmic com- 
plexity and challenge the current model of program devel- 
opment in this research field. In contrast to other research 
areas, in which large "community codes" are being used, 
the field of strongly correlated systems has so far been based 
mostly on single codes developed by individual researchers 
for particular projects. While the simple algorithms used 
a decade ago could be easily programmed by a beginning 
graduate student in a matter of weeks, it now takes sub- 
stantially longer to master and implement the new algo- 
rithms. Thus, their use has increasingly become restricted 
to a small number of experts. 

The ALPS project aims to overcome the problems posed 
by the growing complexity of algorithms and the special- 
ization of researchers onto single algorithms through an 
open-source software development initiative. In release 1.3 
new features desired by ALPS users, especially experimen- 
talists, have been added, as well as a DMRG program for 
interacting problems. 

To quote from the previous publication [1], the goals of 
the ALPS project are to provide: 

- standardized file formats to simplify exchange, dis- 
tribution and archiving of simulation results and to 
achieve interoperability between codes. 

generic and optimized libraries for common aspects 
of simulations of quantum and classical lattice models, 
to simplify code development. 

- a set of applications covering the major algorithms. 

- license conditions that encourage researchers to con- 
tribute to the ALPS project by gaining scientific credit 
for use of their work. 

- outreach through a web page [2], mailing lists and 
workshops to distribute the results and to educate re- 
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searchers both about the algorithms and the use of the 
applications. 

- improved reproducibility of numerical results by 
publishing source codes used to obtain published results. 

- an archive for simulation results. 

The ready-to-use applications are useful both for theo- 
reticians who want to test theoretical ideas about quantum 
lattice models and to explore their properties, as well as 
for experimentalists trying to fit experimental data to theo- 
retical models to obtain information about the microscopic 
properties of materials. 

In the following, we present a quick review of these as- 
pects of the ALPS project, focusing on new features in re- 
lease 1.3. 

2. File formats 

The most basic part of the ALPS project is the definition 
of common standardized file formats suitable for a wide 
range of applications. Standardized file formats enable the 
exchange of data between applications, allow the develop- 
ment of common evaluation tools, simplify the application 
of more than one algorithm to a given model, and are a 
prerequisite for the storage of simulation data in a com- 
mon archive. The ISO standard XML [13] was chosen for 
the specification of these formats because it has become the 
main text-based data format on the internet and because 
it is supported by a large and growing number of tools. 

We have designed a number of XML schemas [14] to 
describe 

the input of simulation parameters, 

- the lattices, 

- quantum lattice models, 

- and the output of results. 

A detailed specification of the formats is given on our web 
pages [2,14]. As examples we show a lattice and a model 
definition in figures 1 and 2. Any of the ALPS applications 
can be run by providing an input file, specifying the sim- 
ulation parameters, together with lattice and model defi- 
nitions (figures 1 and 2) to that application (provided the 
application supports that type of model) . Standardized for- 
mats that extend across all applications reduce the learning 
curve for using the applications and allow common tools to 
be used to analyze the data. 

2.1. New features 

2.1.1. Lattice definitions 

New features in the lattice definitions in release 1.3 in- 
clude the possibility to specify inhomogeneity and deple- 
tion. While in a regular lattice, the model Hamiltonian is 
the same for all bonds or sites with the same type, in an 
inhomogeneous lattice the couplings in the model can be 
different for each bond or site. Examples are disordered sys- 
tems, with randomly chosen couplings such as spin glasses, 
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<LATTICEGRAPH name = "depleted inhomogeneous square lattice"> 
<FINITELATTICE> 

<LATTICE dimension="2"/> 

< EXTENT dimension="l" size="L"/> 

< EXTENT dimension="2" size="L"/> 

<B0UNDARY type="periodic"/> 
</FINITELATTICE> 
<UNITCELL> 

<VERTEX/> 

<EDGE> 

<S0URCE vertex="l" off set="0 0"/> 

<TARGET vertex="l" off set="0 l"/> 
</EDGE> 
<EDGE> 

<S0URCE vertex="l" off set="0 0"/> 
<TARGET vertex="l" offset="l 0"/> 
</EDGE> 
</UNITCELL> 

<INHOMOGENEOUSXVERTEX/x/INHOMOGENEOUS>} 
<DEPLETION> 

< VERTEX probability="DEPLETION" seed="DEPLETION_SEED"/> 
</DEPLETION> 
</LATTICEGRAPH> 

Fig. 1. The definition of a square lattice with one site (vertex) per 
unit cell and bonds (edges) only to nearest neighbors. First the 
dimension, extent, and boundary conditions of the Bravais lattice 
arc described in the <FINITELATTICE> element, then the unit cell 
including the bonds in the lattice is defined. New features in release 
1.3. include inhomogeneity and depletion. 

or systems in spatially varying trapping potentials, such as 
optical lattices in harmonic traps. Inhomogcncities can be 
specified for vertices (sites) and edges (bonds), either for 
all (such as for all vertices in figure 1), or only for one type 
of vertex or edge. 

The other new feature is depiction of a lattice, where 
a fraction of sites or bonds is randomly removed from a 
lattice. Currently only site depletion is implemented, but 
additional types of depletion can be added easily if the 
need arises [15]. In the example in figure 1, the fraction of 
depleted sites is specified to be passed in the DEPLETION 
input parameter, and the random number generator seed 
DEPLETTON_SEED can be changed to give different random 
realizations of the depletion pattern. 

2.1.2. Model definitions 

Figure 2 shows the definition of the Hamiltonian of a 
bosonic Hubbard model in a harmonic trap 

i i 
(id) 1 

The first new feature of release 1.3, shown in this example, 
is the ability to specify composite site and bond operators, 
such as the double_occupancy site term or the bosonJiop 
bond term. These site and bond operators can be used in 
the specification of the model and of measurements. 

The other new feature, also shown in the example, is the 
ability, in conjunction with the specification of an inhomo- 



<BASIS name="boson"> 
<SITEBASIS> 

<PARAMETER name="Nmax" def ault="inf inity"/> 

<QUANTUMNUMBER name="N" min="0" max="Nmax"/> 

<0PERAT0R name="bdag" matrixelement="sqrt (N+l) "> 
< CHANGE quantumnumber="N" change="l"/> 

</0PERAT0R> 

<0PERAT0R name="b" matrlxelement="sqrt (N) "> 

< CHANGE quantumnumber="N" change="-l"/> 
</0PERAT0R> 

<0PERAT0R name="n" matrlxelement="N"/> 
</SITEBASIS> 
</BASIS> 

<SITE0PERAT0R name="double_occupancy" site="x"> 

n(x)*(n(x)-l)/2 
</SITE0PERAT0R> 

<B0ND0PERAT0R name="boson_hop" source="x" target="y"> 

bdag (x) *b (y) +bdag (y ) *b (x) 
</B0ND0PERAT0R> 

<HAMILTONIAN name="trapped boson Hubbard"> 
<PARAMETER name="mu" def ault="0"/> 
PARAMETER name="t" def ault="l"/> 
PARAMETER name="U" def ault="0"/> 
<PARAMETER name="K" def ault="0"/> 
<BASIS ref="boson"/> 
<SITETERM slte="i"> 

-mu*n(i) +U*double_occupancy(i) + K/2* ( (x-L/2) ~2 + (y-L/2)~2) 
</SITETERM> 

<BDNDTERM source="i" target="j"> 

-t*boson_hop(i , j ) 
</B0NDTERM> 
</HAMILTONIAN> 

Fig. 2. The definition of a bosonic Hubbard Hamiltonian (1) in a 
harmonic trap. After describing the local basis for each site and 
operators acting on it, the first new feature in release 1.3 is the ability 
to define composite site and bond operators, such as the double 
occupancy and bosonic hopping terms, which are then used in the 
definition of the Hamiltonian. 

geneous lattice, to use site-dependent couplings, such as 
the harmonic trapping potential K which depends on the x 
and y coordinates of the site. This feature has been used for 
the simulation of cold bosonic gases in optical lattices [16]. 

2. 2. Future plans 

The immediate goal for release 1.4 will be to add multi- 
site interaction terms and measurements, such as ring ex- 
change terms. 

With the experience gained in using common file formats 
for several years across a number of different applications, 
we next plan to design more widely used common formats, 
such as standard formats for results of Monte Carlo simu- 
lations that have been defined at a workshop in September 
2006. 

3. Libraries 

The ALPS libraries are the foundation of all the ALPS 
applications, providing functionality common to all of 
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them, such as file I/O, expression evaluation, parallel job 
scheduling, or the statistical analysis of Monte Carlo mea- 
surements. These libraries make full use of object-oriented 
and generic programming techniques [18], which allows 
them to be very flexible without losing any performance 
compared to FORTRAN programs. 

3.1. New features 

The new library features are mainly hidden from the 
application users, providing support for the new features 
in the input files, discussed above, or providing support 
for new features in the applications, such as translational 
symmetries in the diagonalization codes, and the ability 
to define custom measurements, which will be discussed 
below. 

One new feature that should be mentioned explic- 
itly is the inclusion of random number generators in the 
expression library. Expressions for couplings and other 
parameters can make use of new functions randomO, 
normal_random ( ) , and integer_random(n) to obtain 
uniform random numbers in the interval [0, 1), normally 
distributed random numbers or integer random numbers 
between and n. This feature is most useful in conjunction 
with inhomogeneous lattices to define disordered models. 

Another change to be mentioned is that since the pars- 
ing of command line options to the application is now han- 
dled by the Boost program options library [17], multi-letter 
command line options now need to be preceded by two hy- 
phens as in — Tmin, instead of the previous usage -Tmin. 

3. 2. Future plans 

Code from the ALPS libraries has recently been used in 
the development of Boost [17] libraries. These include op- 
timizations to the Boost. Serialization library, new features 
for the Boost. Random library, a proposed Boost. MPI li- 
brary and a Boost. Accumulator library for statistical esti- 
mates. Once accepted in the peer review process, and pub- 
lished in a future Boost release, the ALPS codes will be 
changed to use these Boost libraries. The ALPS libraries 
will then contain less general purpose codes (which will have 
been moved to Boost) and will be more focused on domain- 
specific libraries for the simulation of quantum models. 

4. Applications 

In addition to common libraries, the ALPS project in- 
cludes a number of ready to use applications implement- 
ing the most important unbiased algorithms for quantum 
many body systems. The applications all share the same 
file formats, simplifying their use, reducing the learning 
curve, and enabling the easy investigation of a model with 
more than one method. Tutorials on the use of the appli- 
cations are included with the sources that can be found on 
the ALPS web page [2]. 



4-1- Exact diagonalization 

The existing exact diagonalization [19,11] programs 
sparsediag and fulldiag have been substantially im- 
proved, sparsediag calculates the ground state and low 
lying excited states of quantum lattice models using the 
Lanczos [20] algorithm, while fulldiag calculates the 
complete spectrum of quantum lattice models and from it 
all thermodynamic properties. 

The first major improvement is the use of translation 
symmetry to reduce the Hilbert space dimension. This 
speeds up the calculations and allows larger systems to be 
calculated. In addition, by calculating the energy eigen- 
values separately for each momentum, the momentum- 
resolved excitation spectrum can be calculated. 

The second new feature is the possibility to define cus- 
tom measurements of averaged or local site- and bond op- 
erators, as well as arbitrary 2-point correlation functions. 
As an example, for the inhomogeneous bosonic Hubbard 
model of figure 2, one can specify the measurement of the 
average double occupancy, local density, density correla- 
tions, and Green's function by defining the following input 
parameters, which use the definitions of site and bond op- 
erators provided with the model: 
MEASURE_AVERAGE [Double] = double_occupancy 
MEASURE_LDCAL [Local density] = n 
MEASURE_CDRRELATIDN [Density correlation] = n 
MEASURE_CDRRELATIDN [Green function] = "bdag:b" 
MEASURE_STRUCTURE_FACTDR [Spin Str Fact] = Sz 

In the future we plan to support point group symmetries 
and to add the calculation of dynamical correlation func- 
tions. 

4-2. Classical Monte Carlo 

The Monte Carlo program for classical magnets employ- 
ing local and cluster updates [21] has been extended ac- 
cording to wishes of several users. The most notable new 
features are the possibility to specify single ion anisotropies 
and anisotropic interactions. 

4-3. Quantum Monte Carlo 

The quantum Monte Carlo programs have also been up- 
dated in release 1.3. 

The "looper" program using the loop cluster algorithm 
[7,22] in a stochastic series expansion [23] (SSE) and path- 
integral representation for quantum magnets has been re- 
placed with a new version, offering 40-80% speed increase 
and supporting longitudinal and transverse magnetic fields 
as well as single ion anisotropies. 

The directed loop QMC program [24] in an SSE repre- 
sentation for bosonic and magnetic quantum lattice models 
has been extended to offer also the "Locally optimal worm 
algorithm" (LOWA) [25] updates. 

The worm i ilgorithm [26] program has been updated to 
allow nonlocal interactions [27] . 
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All QMC program have been extended by new measure- 
ments, including correlation functions, structure factors, 
and stiffness or superfluid density. In addition, the same 
custom measurements as discussed above for exact diago- 
nalization can be used in these codes, as long as the mea- 
surements are diagonal in the chosen basis. The looper pro- 
gram additionally measures local susceptibilities and trans- 
verse magnetization. 

4-4- Density Matrix Renormalization Group 

In addition to the DMRG [9,10,11] program for nonin- 
teracting particles (an adaptation of the program in Refs. 
[28,29] to the ALPS libraries), release 1.3 contains a new 
DMRG program for interacting many body systems, with 
the ability to calculate local quantities and two-point cor- 
relation functions. 

4-5. Tools 

Release 1.3 adds a database tool, to archive a collection of 
simulation results, and to allow fast and efficient extraction 
and evaluation of results. The tool, based on the SQlite [30] 
database is easy to use, based on a simple flat file database 
and does not require to set up a complex database appli- 
cation. 

5. License 

The ALPS software is distributed under the ALPS li- 
brary [31] and application [32] licenses, which have re- 
mained unchanged. The use of the codes is free for non- 
commercial applications but carries a citation requirement 
to publications, such as the current one, describing the li- 
braries and codes. In addition, besides the required citation 
of the publication describing the ALPS application, users 
are encouraged to cite the publication describing the algo- 
rithm used in the applications. 

6. Distribution, outreach development 

The ALPS web page [2] is the central information reposi- 
tory of the ALPS collaboration. It makes source codes of the 
ALPS libraries and applications available, provides access 
to the documentation, and distributes information about 
the mailing lists, new releases, and workshops. 

For this release we have replaced the static web page 
by a Wiki system [33] , which enables developers and users 
to easily and quickly add new contents and improve the 
instructions and tutorials on the web page. 

Another new feature is cdALPS [34] by G. Pawlowski, a 
bootable Knoppix Linux ALPS CD, which allows to run the 
ALPS tutorials and simple ALPS codes on a PC without 
installing the full ALPS software. 

Development of the ALPS libraries and applications is 
coordinated through both mailing lists and semi-annual 
developer workshops. Interested researchers are invited to 



join the mailing list and participate in the workshops. 
7. Future plans 

7.1. New application 

The following new applications are planned for future 
releases: 

- A series expansion [35] code for both perturbation 
series and high temperature series in release 1.4 

- Also in release 1.4, the addition of generalized sam- 
pling algorithms, such as the quantum version [8] of 
the Wang-Landau algorithm [36] and optimized en- 
semble methods [37] to the quantum Monte Carlo 
algorithms to obtain thermodynamic quantities over 
large temperature ranges. 

- Implementation of parallel tempering [38] for the clas- 
sical and quantum Monte Carlo algorithms including 
feedback-optimized temperature sets [39] . 

- A QMC algorithm in the valence bond representation 
for quantum spin systems [40] 

- A continuum QMC program for bosons, implement- 
ing the continuous space worm algorithm [41]. 

7. 2. Realistic model: interface to band structure codes 

Also planned for release 1.4 is an interface of the ALPS 
applications to band structure codes. Instead of running 
simulations on phenomenological toy models, such as a sim- 
ple Hciscnbcrg model with only a few different couplings, 
it will be possible to run the ALPS applications on a com- 
plex, realistic model obtained from ab-initio band struc- 
ture calculations. In the past this was done manually, e.g. 
in reference [42], the exchange couplings determined in an 
LDA+U calculation [43] using the Lichtenstein method [44] 
were used to manually prepare input files for QMC simu- 
lations using the loop algorithm [7] . 

Release 1.4 of ALPS will contain tools to automate such 
calculations, based on a standard XML I/O format for 
band structure codes developed by T. Schulthess and M. 
Summers at the Oak Ridge National Laboratory. As an ex- 
ample of such an automated calculation using a prototype 
tool by A. Kozhevnikov, we show in figure 3 a compari- 
son of experimental measurements on the spin ladder com- 
pound SrCu 2 03 [45] with automated ab-initio simulations 
("LDA+U+ALPS"). 

7.3. Dynamical mean field theory (DMFT) solvers 

In a further stage, we plan to expand the ALPS project to 
include dynamical mean field theory (DMFT) [46] solvers. 
We have designed a DMFT framework for QMC solvers 
and implemented the standard Hirsch-Fye solver [47], as 
well as two new continuous time solvers [48,12]. The first 
application of the framework, a careful comparison of the 
performance of the three solvers [47,48,12] is in preparation 
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Fig. 3. Comparison of experimental measurements on the spin lad- 
der compound SrCu2 03 [45] with automated ab-initio simulations 
using "LDA+U+ALPS". The deviation at low temperature is due 
to impurities in the experimental sample, causing a Curie-Weiss 
like divergence at low temperatures. The dominant couplings are 
an intra-chain coupling J = 162meV and an inter-chain coupling 
J ± = 76meV 

[49]. Work is in progress to expand the framework to in- 
clude ED and DMRG based solvers, and we expect a public 
open-source release of the first codes soon. Interested re- 
searchers can contact us to participate in the development 
or to obtain pre-release versions of the codes. 

8. Conclusions 

The ALPS project is continuously evolving. Researchers 
interested in announcements of new releases, information 
about workshops, or in contributing to the ALPS project 
are encouraged to sign up to the mailing lists on our web 
page [2]. 

We acknowledge support by the Swiss National Science 
Foundation, the Kavli Institute for Theoretical Physics at 
the University of California at Santa Barbara, the Aspen 
Center for Physics, CNPq (Brazil) , the Network for Interna- 
tional Development and Cooperation (NIDECO) at ETH 
Zurich, the RGC grant of Hong Kong, and the Deutsche 
Forschungsgemeinschaft through SFB 602. 
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